For Learning Purpose Only 

%Question 1 part a 

xmi n = -20; 
xmax = 30; 
ymi n = ; 
ymax = 8 ; 
n=- 20:30 ; 
t = 0:30; 

% Creating 31 impulses from n=0 (origin) to n=30 
IMPulses = eye (31) ; 

% Initializing variable for Step Function. 
STEp = zeros (1,31) ; 

% (Now) Adding Impulses to create a step. 

for i=l:31 

STEp = STEp + IMPulses (i, :) ; 

End 

s tern (t, STEp) ; 

axis([xmin xmax ymin ymax]) 



Do not submit unless you understand the code you are submitting. 



For Learning Purpose Only 

% Question No. 1 part b 

xmi n = -20; 
xmax = 30; 
ymi n = ; 
ymax = 10; 
n=- 20:30 , ■ 
t=0:30; 

LENgth = input ('Enter Desired Length : ') ; 
x = round (10* rand ( 1, LENg th )) ; 

disp( 'Generated Impulses are : '); 

disp(x) ; 

display ( 'Press any key to Continue'); 

pause 

IMPulses = eye (31) ; 
SIGnal = zeros (31); 
DESiredSignal = zeros(l,31); 

for i=l : LENgth 

SIGnal (i,:) = x (i) *IMPulses (i , : ) ; 

end 

for j=2:31 

SIGnal (1,:) = SIGnal (1,:) + SIGnal (j ,:) ; 

end 

DESiredSignal = SIGnal (1,:) ; 

s tern (t, DESiredSignal) 
axis([xmin xmax ymin ymax]) 



Do not submit unless you understand the code you are submitting. 



For Learning Purpose Only 

function OUTput = CONvolutionl (x,y) 



LENgth = length (x) + length (y) - 1 ; 
DS = zeros (1, LENgth) ; 
eDS = DS; 

if (length (x) >length (y) ) 

Osmall = y; 

Olarge = x; 

else 

Osmall = x; 

Olarge = y; 

end 

for j = 1 : length (Osmall) 

DS(j : (length (Olarge) +j-l) ) = Osmall (j ) *01arge ; ■ 

eDS (j : (length (Olarge) +j-l) ) = eDS (j : (length (Olarge) +j -1) ) + 

DS(j : (length (Olarge) +j-l) ) ; 
end 

OUTput = eDS; 



Do not submit unless you understand the code you are submitting. 



For Learning Purpose Only 

%Question 2 part b 

hi = 4*[0 0.5 0.25 0] / 

h2 = [1 -4 0] ; 

h3= [0000010]/ 

h4 = h3; 

h5 = h2; 

h34 = CONvolutionl (h3,h4) ; 

hl34 = CONvolutionl (hi, h34) ; 

hl2 = CONvolutionl (hi, h2) ; 

%##### this could be a function ######### 

if (length (hl34) >length (hl2) ) 

LENgth = length (hl34) ; 

PEL1 = (-1) * hl34; 

PEL2 = hl2; 

PEL1 (1 : length (PEL2) ) = PEL1 (1 : length (PEL2) ) + PEL2; 

else 

LENgth = length (hi 2 ) ; 

PEL1 = hl2; 

PEL2 = (-1) * hl34; 

PEL1 (1 : length (PEL2) ) = PEL1 (1 : length (PEL2) ) + PEL2; 

end 

hl234=PELl; 

%##### this could be a function ######### 

if (length (hl234) >length (h5) ) 

PEL3 = hl234; 

PEL4 = h5; 

PEL3(1 : length (PEL2) ) = PEL1 (1 : length (PEL2) ) + PEL2; 

else 

PEL3 = h5; 

PEL4 = hl234; 

PEL3(1 : length (PEL2) ) = PEL1 (1 : length (PEL2) ) + PEL2; 

end 

hl2345=PEL3; 

t = 0: (length (hl2345)-l) ; 

stem(t,hl2345) / 



Do not submit unless you understand the code you are submitting. 



For Learning Purpose Only 

% Question No. 3 

phasesl = [ (0) (pi/4) (pi/2) (3*pi/4) (pi) ] ; 
ampsl = [0.5 0.6 0.7 0.8 0.9]; 
freqsl = [0.02 0.01 0.4 0.2 0.005]/ 
n = 0: .1:100; 

X = zeros (1 , length (n) ) ; 

for i=l : 5 

X = X + ampsl (l,i) * (sin( (freqs(l,i) *n) + phasesl (l,i) ) ) ; 

end 

t = 0: (length (X)-l) ; 

plot(n,X) ; 



Do not submit unless you understand the code you are submitting. 



